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Dynamics of Wilson loops in pure Yang-Mills theories is analyzed in terms of random walks of 
the holonomies of the gauge field on the gauge group manifold. It is shown that such random 
walks should necessarily be free. The distribution of steps of these random walks is related to the 
spectrum of string tensions of the theory and to certain cumulants of Yang-Mills curvature tensor. 
It turns out that when colour charges are completely screened, the holonomies of the gauge field 
can change only by the elements of the group center, which indicates that in the screening regime 
confinement persists due to thin center vortices. Thick center vortices are also considered and the 
emergence of such stepwise changes in the limits of infinitely thin vortices and infinitely large loops 
is demonstrated. 
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I. INTRODUCTION 

Wilson loops are the most popular order parameters 
for Yang-Mills theories without dynamical quarks. Wil- 
son loop Wr [C] is defined as the expectation value of 
the trace of the holonomy of the gauge field over the 

loop C: Wr [C] = ( Tr rV exp (^i J dx^Af,^ ) 1]. Inter- 
action potential Vr (r) between two static colour charges 
which transform under some irreducible representation 
R of the gauge group can be measured using the Wil- 
son loop Vl/i{[Crxt], where Crxt is a rectangular loop 
of size r X t with t co. As the Wilson loop is 
in fact the amplitude of propagation of static colour 
charges along the loop C [3], it is related to the poten- 
tial Vr (r) as Wr [Crxt] = exp {—Vr (r) t). It is known 
that in the confining phase of the theory static colour 
charges are connected by a thick chromoelectric string 
with constant tension aR, which gives rise to linear in- 
teraction potential Vr (r) = ctr r at sufficiently large 
distances. Correspondingly, for sufficiently large loops 
Wilson loops Wr [C] decay exponentially with the mini- 
mal area S [C] of the surface spanned on the loop C, i.e. 
Wr [C] = exp {—aR S [C]). Contribution of charge self- 
energies shows up in the dependence of Wr [C] on the 
perimeter of the loop C. 

Thus confining properties of the theory are character- 
ized by the whole spectrum of string tensions cjr. The de- 
pendence of aR on the representation of the gauge group 
gives some additional information on the properties of 
Yang-Mills theory. For instance, lattice simulations show 
that at intermediate distances (0.2 — Ifm for SU{3) gauge 
group) aR is proportional to the eigenvalue C'2 r of the 



second-order Casimir operator in the representation R 
(Casimir scaling) 0, 01 . Casimir scaling is naturally ex- 
plained in the framework of stochastic vacuum models 
[j]. At very large distances Casimir scaling should be 
violated because of screening of colour charges by gluons 
Q . The reason is that when colour charges are separated 
by sufficiently large distance, it costs less energy to create 
a bound state of a colour charge and some number of glu- 
ons than to create confining string between bare colour 
charges. For instance, when SU{N) colour charges are 
completely screened, string tension depends only on the 
iV-ality of the representation 0, 0| and is proportional 
to the lowest eigenvalue of the second-order Casimir oper- 
ator among all irreducible representations with the same 
iV-ality. 

An interesting description of the spectrum of string 
tensions of Yang-Mills theories was proposed recently 
in [1, where all information about Wilson loops 

Wr [C] in irreducible representations of the gauge group 
was encoded in a single function on the gauge group man- 
ifold, namely, the probability distribution p[g;C] of the 
holonomies of the gauge field over the loop C. The holon- 
omy g [C] = Vexp J dx^^Afj^ is defined only modulo 

gauge transformations g [C] —> hg [C] h~^, which implies 
that p [g; C] should be invariant under such transforma- 
tions, i.e. p [g;C] ^j}[hgh-^;C]. Thus p[g;C] should 
be defined as [8j,l9j,l3: 



p[g-C] = {S, {g,g[C])) 



(1) 



where Sc{g,g') = J2xR{g)XB.{g') is the delta-function 
R 

on the group classes and XR {9) are the group characters. 
Using the character expansion of 5c {g,g'), one can im- 
mediately express p[g;C] in terms of the Wilson loops 
Wr[C]: 
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p[g;C]^Y.^'^(9)WR[C] 
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satisfies the free diffusion equation: 
d 



dS [C] 



p [g; C] = aoAp [g; C] 



(4) 
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FIG. 1: Small increment of the area of the minimal surface 
spanned on the loop C 



Conversely, all Wilson loops can be calculated if the 
probability distribution p[g;C] is known: Wr [C] ~ 
J dg XR (5) P b; C*]. ft can be shown that p [g; C] is the 
Wilson loop in the regular representation of the gauge 
group 

An advantage of such description is that the evolution 
of p [g; C] can be interpreted in terms of a random walk 
of the holonomy g [C] on the gauge group manifold, the 
properties of this random walk being directly related to 
the spatial distribution of non-Abelian magnetic flux, fn 
a particular configuration of gauge fields, the holonomy 
g [C] moves along some path on the group manifold as the 
area of the loop is gradually increased. When one cal- 
culates expectation values of Wilson loops in quantum 
theory and sums over all field configurations, weighted 
sum over all such paths which start in the group identity 
and end in the element g yields the probability distri- 
bution p[g;C]. In the case of pure Yang-Mills theory 
in Euclidean space-time all field configurations and, con- 
sequently, all such paths are summed with nonnegative 
weights, therefore such weighted sums over paths on the 
group manifold can be described as random walks on Lie 
groups The position of random walker corre- 

sponds then to the holonomy g [C] modulo gauge trans- 
formations g [C] — > hg [C] h~^. If the loop C is slightly 
deformed in such a way that the area of the minimal sur- 
face increases from S [C] to S [C] = S [C] + SS (see Fig. 
[T|), the holonomy g [C] changes as: 

g [C] = g [C] (I + ih (xo, x) F^, (x) h (x, xq) SS^ (3) 

where h{x,y) = "Pexp J dx^Af^^ is the multiplica- 
tive integral of the gauge field over the segment of the 
loop C bounded by the points x and y. The point 
Xo is the initial point used to calculate the holonomy 

g[C]: g[C] =7'exp[z / dx'^A^]. The combination 

\ C;xo J 

F^u{x) — h{xQ,x) F^^ {x) h^^ {x,xq) is usually called 
the shifted curvature tensor Q. Thus each step of the 
random walk of g \C] corresponds to a small change of 
non-Abelian magnetic flux through the loop C . 

One of the basic observations of f8, 9, is that when 
Casimir scaling holds, the probability distribution p[g;C] 



where A is the Laplace operator on the group man- 
ifold and Co is some constant. Random walk which 
is described by ^ is simply the Brownian motion on 
the group manifold, which consists of a large number 
of small statistically independent random steps. Corre- 
spondingly, the distribution of non-Abelian flux through 
the loops of intermediate sizes is almost random, which 
fits nicely in the model of stochastic vacuum 

Free diffusion equation ([4]) and the corresponding ran- 
dom walk should be somehow modified in order to re- 
produce screening effects. In [1, Q random walk in Zff- 
symmetric potential was considered and it was shown 
that diffusion in such potential reproduces transition 
from Casimir scaling to complete screening due to mix- 
ing of different representations with equal A'^-alities. An 
advantage of such approach is that one can solve the dif- 
fusion equation in terms of path integral and obtain the 
effective action for Wilson loops. However, this approach 
is to a large extent phenomcnological - in particular, the 
choice of the potential which breaks SU [N) symmetry 
of the equation down to Z^v is to a large extent arbitrary 
and affects only the transition from Casimir scaling to 
scree ning . 

In the most general form of the diffusion equation 
for p [g] C] was obtained using the classical loop equations 
and the cumulant expansion theorem jlSi] . It turned 
out that this diffusion equation should have the form of 
the so-called Kramers-Moyall cumulant expansion, which 
describes a free random walk with an arbitrary distri- 
bution of steps. Free random walk is understood here 
as a random walk with position-independent step distri- 
bution. Thus exact equation for p [g; C] cannot include 
potential or drift terms. 

The aim of this paper is to consider the most general 
free random walk on the group manifold and to relate 
its properties to the spectrum of string tensions of Yang- 
Mills theory. Although in general it is not possible to find 
the effective action which describes a free random walk, 
the distribution of steps of the random walk of g [C] has 
a more natural physical interpretation than an external 
potential in the diffusion equation. In particular, it will 
be shown that the distribution of steps of the random 
walk changes dramatically when Casimir scaling changes 
to screening. While the free diffusion equation @ reflects 
stochasticity of Yang-Mills vacuum at small and interme- 
diate distances, random walk which describes diffusion of 
Wilson loops in the screening regime can only consist of 
discrete jumps by the elements of the group center, which 
can be naturally explained in terms of thin center vor- 
tices [H, [13, [2l| . It will be also demonstrated that in the 
case of center vortices with finite thickness the holonomy 
g [C] can change by any group element, but the distri- 
bution of steps is still peaked near the elements of the 
group center, the widths of the peaks being proportional 
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to vortex thickness and inversely proportional to the area 
S [C] of the minimal sm'face spanned on the loop. 

The structure of the paper is the following: in the sec- 
tion [IT] a general case of a free random walk of g [C] on 
the group manifold is considered. The distribution of 
steps of such random walk is expressed in terms of the 
spectrum of string tensions of the theory and in terms 
of the cumulants of the shifted curvature tensor i^^^. In 
the section Hill the random walk of g [C] in the screening 
regime is investigated and it is shown that in the case 
of complete screening the holonomy g [C] is only allowed 
to change by the elements of the group center. In the 
section IIVI such discrete jumps are interpreted in terms 
of thin center vortices. Vortices of finite thickness are 
considered in the section |Vl where it is shown how such 
jumps by the elements of the group center arise in the 
limits of infinitely thin vortices or infinitely large loops. 

II. FREE RANDOM WALKS ON THE GROUP 
MANIFOLD 

In [ll| it was shown that the most general form of the 
diffusion equation for the probability distribution of the 



holonomies g [C] on the gauge group manifold is: 



— p[g;C]=Y^ ir-^- [C] Va, . . . Va, P [ff; C] (5) 



where 7^°i-- °fe [(7] are some coefficients which can in gen- 
eral depend on the loop C and Vq are the generators of 
left shifts on the group manifold. The action of Vq is de- 
fined by the identity (1 + e'^Va) f (g) = f ((1 + le^Ta) g), 
where are arbitrary infinitely small parameters, Ta are 
the generators of the gauge group and / (g) is an arbi- 
trary function on the group manifold. The term with 
first-order derivative is prohibited in ([5]) by gauge invari- 
ance. It was shown in that the coefficients t]"'^ -"-'' [C] 
are directly related to the cumulants of the shifted cur- 
vature tensor F^^, which are the basic objects in the 
method of field correlators [3] ■ This relation can be for- 
mally written as (ll| : 



SIC] S[C] 

where the integration is performed over the surface of the minimal area S [C] spanned on the loop C. It should 
be noted that for any point x the shifted curvature tensor f^^, (x) transforms under gauge transformations as the 
curvature tensor in some fixed point xq, therefore the cumulants of i^^"^^ (x) are well-defined when all xi, . . . ,Xk are 
different. Gauge invariance implies only that the cumulants of F^^^ (x) and the coefhcients j^^i [C*] should be colour 
singlets. In particular, all ri"'^---"-^ ^C] with odd k should vanish. 

The equation ([5]) is a general equation which describes an arbitrary free random walk on the group manifold, i.e. 
a random walk with position-independent distribution of steps. If P [g'] C] SS is the probability that g [C] changes 
to g [C] — g'g [C] as the area of the surface spanned on the loop C increases from S [C] to S [C] ~ S [C] + 6S, the 
evolution of the probability distribution p [g; C] is described by the following equation: 

p [g; C]^ J dg'P [g'g-';C] p [g'; C] - p [g; C] J dg'P [g'; C] = 

= J dg' (^P [g'g~';C] ~ S {g', g) j dg" P [5"; C]) p [g'- C] (7) 

where 5{g',g) is the delta-function on the group manifold. Although the distribution of steps in ([7]) can in general 
depend on the loop C, this dependence can be neglected if the Wilson area law holds exactly and if the spectrum of 
string tensions is known. Using the definition ([T|) the following integral equation can be obtained for p [5; C]: 

p[g-C]^~Y. ^« <yRWR [C] - - ^ XR (5) <y (R) [ dg'p W; C] XR (5') = 
= -JdhJdg'J2 dRXR {ghg'-'h-') a {R)p[g'; C] ^ J dg' U ^ dnXR {gg'-') <r (R)\ p [g'; C] (8) 

R \ R / 



where the identity / dhxR {ghfh ^) — dp^XR (ff) X-R (/) tions g [C] — > hg [C] h ^ were used. Comparing the equa- 
and the invariance ol p[g;C] w.r.t. gauge transforma- 
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tions ([8]) and ([7]), one can finally express the step distri- 
bution P [g;C] in terms of string tensions an: 

P[5;q = -^dflaflXfl(5) (9) 

R 

If string tensions depend on the distance between colour 
charges, for instance, because of screening, P [g; C] 
should also depend on the loop C. This dependence will 
be investigated in the section |V] using the model of thick 
center vortices. It is also interesting to note that the pos- 
itivity of P [g; C] for g 7^ 1 is some constraint on possible 
spectrum of string tensions of the theory, which follows 
from the positivity of path integral weight for pure Yang- 
Mills theory in Euclidean space-time. 

The equation ([5|) is simply the gradient expansion of 
the general equation ([7]). Indeed, p{g') in (O can be re- 
placed by p [g'g-'g; C] = exp {-^ {99'-') S7a)p [9; C], 
where the function ^" (g) is defined by the following iden- 
tity: 

e^piiCi9)Ta)=g (10) 

Expanding the exponential in powers of {g9'~^), one 
can express the coefficients 7y"i [C*] in terms of the 
step distribution P[g;C]: 

V"'-^" [C] = ^ j dgP [g- C] (g) . . . ^ ^ (5) (11) 

As all 7y°i- -'''= [C] are colour singlets, the distribution of 
steps should satisfy P [hgh~^; C] = P [g; C]. 

For example, if Casimir scaling holds, (Tr — <tqC2 r, 
P [g; C] = cToAS {g, 1) and only the coefficient r;"'' [C] = 
aoSab is not equal to zero. Correspondingly, only the 
second-order cumulant of the shifted curvature tensor 
is nonzero, which corresponds to the limit of gaussian- 
dominated stochastic vacuum Nl. 



III. COMPLETE SCREENING AND DISCRETE 
JUMPS ON THE GROUP MANIFOLD 

In the case of complete screening string tensions 
depend only on the A^-ality of the representation R: 
<^R = <^ {^r)- A^-ality Ar of the representation R charac- 
terizes its transformation properties w.r.t. the elements 
of the group center Zjq = {e"^ '^j, fc = 0, 1, . . . , iV — 1: 
XRigz) = XRig)^^"'^ ^ £ Zjsi. A simple calculation 
shows that if gr — a (A/j), P [g; C] is given by a finite 
sum over the elements of the group center: 

P[9\C]^ E ^i^)H9.z) (12) 

where (z) is minus the Fourier transform of the spec- 
trum of string tensions w.r.t. the group center: 

N-l 

77(z) = -iV-i5]a(A)zA (13) 

A=0 



Thus when colour charges are completely screened, the 
evolution of the probability distribution p[g;C] is de- 
scribed by the following equation: 

^^^p[g-C]=Y.Tl{z)p[zg-C] (14) 

This equation implies that the random walk of g [C] 
consists of jumps by the elements of the group center 
only. 77 (z) with z ^ 1 is the probability of a jump by z 
per unit area. As the string tension between charges with 
zero A^-ality should vanish due to complete screening by 
gluons 0, ^ ''7 (-2) = 0, which is simply the conserva- 

tion of probability flow in the language of random walks. 

Thus the distribution of steps of the random walk of 
g [C] is completely different in Casimir scaling and in 
screening regimes. When Casimir scaling holds, random 
walk of g [C] consists of a large number of statistically in- 
dependent small steps, which indicates that non-Abelian 
magnetic flux which penetrates the loop C fluctuates ran- 
domly. On the other hand, the equation ^4]) implies that 
the cumulants in ^ are saturated by some singular field 
configurations. 

IV. THIN CENTER VORTICES 

It is straightforward to guess what is the structure of 
singular field configurations which correspond to discrete 
jumps described by (|14p . By definition the holonomy 
g [C] changes by the element of the group center when 
the loop C is crossed by center vortex [16|. As follows 
from equation g{C\ changes stepwise, which can 

only be explained if the vortices are mathematically thin 
surfaces which are distributed in space with some finite 
density. Correspondingly, 77 (z) 5S is the probability to 
cross thin center vortex which carries magnetic flux z as 
the area of the loop is increased from S [C] to S [C] +5S. 

The statement about the singular field configurations 
can be formulated more precisely using the equations ^ 
and (fTTj) . which relate the step distribution P [g; C], the 
coefficients 77"! ■■■'^fc [C] and the cumulants of the shifted 
curvature tensor of Yang-Mills fields. The coefficients 
^ai...ak YJ^^ calculated if the field configurations 

which contribute to the cumulants in ([6]) are known. On 
the other hand, these coefficients can be found from the 
step distribution P[g]C]. In this section it will be shown 
that the coefficients jy^i - ^fc [C] obtained from the dis- 
tribution (fT2|) agree with the result obtained from ^ if 
center vortices are infinitely thin non-interacting random 
surfaces. This analysis will be extended in the next sec- 
tion, where center vortices of finite thickness will be con- 
sidered and the emergence of the step distribution (fT2|) in 
the limit of infinitely thin vortices will be demonstrated. 

In order to calculate the coefficients rj'^i - '^k [(7] for 
the step distribution (fT2|) . one should calculate the in- 
tegrals of the form / dgS {g, z) (5) ■ • • C'' (s)- Such in- 
tegrals, however, can not be calculated by direct substi- 
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tution g ^ as for the elements of the group center 
the choice of (g) is not unique. In fact, the points 
g — z with z 1 are the only points on the group man- 
ifold where the function (g) is not uniquely defined. 
Indeed, as exp (i^*^ (z) T^) — z, for any group element h 
one has hexp {z)Ta)h-'^ = exp {iS,'' (z) hTah-^) = 
exp (i^^ (z) 0° (h) To) = hzh-'^ = z (here (h) are 
the matrices of the adjoint representation of the gauge 
group). The integrals J dgS {g , z) S^""^ (g) . . . S^"-'' (g) can 
nevertheless be calculated if the delta-functions in (fT2l) 
are defined as the limits of some smooth functions, 
which immediately yields / dgS {g, z) [g) ■ ■ ■ S,""^ (g) = 
/ dh Oil (^) (^) • ■ ■ ^ll (^)- It 



nient to introduce the notation (^""^ (z) . . . (z))^ = 
/ dh Oil (^) ■ ■ • (^) (^)' ^° t^at the coeffi- 

cients Ty^i - '^fc [(7] which correspond to the step distribu- 
tion ([T2|) can be written as: 



The equations ([6]) and (fT5| imply that in this case the 
integrals of the cumulants of the shifted curvature tensor 
should have the following form: 



dS^^^^{x,)...dS^^'''^{xu){{F;i^^{x,)...F;i,^{xu)))= S[C] {az)...eHz))k (16) 



S[C] S[C] 



On the other Jiand, the integrals in (|T6|) can be calculated using the exact expression for the field strength of thin 
center vortex [ISj : 

(^) = r {z) I dE"^ {y) [x - y) (17) 



where S is the vortex surface. The choice of (z) in (1171) is also not unique, and it is usually assumed that physical 
observables should be averaged over all possible (z) jl9l.[20j. It is convenient now to introduce the linking number 
of the vortex surface S and the loop C: 

L[C;S] = e„;3^, j dS^''' {x) J dJ:^^ (y) 5^ {x ~ y) (18) 

SIC] s 

The linking number L [C; E] counts how many times the surface E winds around the loop C. Using the expressions 
(I17p . (jlSp and averaging over all possible vortex configurations and over all (z), one can rewrite the integrals in 
(fT6l) as: 

dS'^^-'^ (xi) . . . dS^^'^" {xk) {{ F^l^^ (xi) . . . F;^^^ (x,) ))=$](( [C; E,] (^) ■ • • (z) )) (19) 

S[C] SIC] "'^^^ 

where E^ denotes the union of the surfaces of all vortices which carry magnetic flux z. As only the cumulants of 
even order are different from zero, all indices ai, . . . , can be pairwise contracted in both (|19p and (flB]) . As for all 
possible (z) the squares (z) (z) are equal, this gives equal constant factors in front of both expressions. After 
that direct comparison of (fT9|) and (fT6|) shows that all even-order cumulants of the linking number L [C; Ez] are equal: 
((L^'^ [C; Ez] )) = 7] (z) [C], ((L^'^+i [C; E^] » = 0. This property is enough to recover the probability distribution of 
the linking number L [C; E^] - it is distributed as the difference of two Poisson-distributed positive integer numbers 
and fc_ with mean values rj (z) S [C] /2 each. The numbers k+ and fc_ can be readily interpreted as the numbers 
of vortices which wind around C leftwards and rightwards. Explicit expression for the probability distribution of 
L [C; Ez] is: 

ML[C;Ez])=exp(-^(z)5[C]) J] (20) 

fc+-fc_=L[C;E^] +' 

I 



The fact that fc_|- and fc_ are distributed by Poisson im- 



plies that thin vortices do not interact, since all intersec- 
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tions of vortex surfaces and the loop C are statistically 
independent events. Recent results of lattice simulations 
indicate that center vortices are indeed mathematically 
thin random surfaces which are distributed in space with 
finite density pvort ~ 24 fm~^ [2l[. A typical vortex 
configuration on the lattice consists of a single perco- 
lating vortex which stretches through the whole physi- 
cal space plus a large number of small vortices, whose 
sizes do not exceed several lattice spacings [2l|, |22| ■ As 
such small vortices can only lead to perimeter-dependent 
effects, the existence of the percolating vortex is cru- 
cial for the area dependence of the expectation values 
{k+) = {k-} =r]{z)S [C] 12. The results obtained in [H 
also give some preliminary indications that the effective 
interaction between thin vortices is rather small. Thus 
the picture of thin non-interacting center vortices may 
be a good approximation for the low-energy dynamics of 
pure Yang-Mills theories. 

V. THICK CENTER VORTICES 

The analysis presented in the previous section is based 
on the assumption that thin center vortices are the only 
field configurations which contribute to the cumulants of 
the shifted curvature tensor in ([6]), which is justified if the 
equation (fT4|) is exact for loops of finite size. However, 
since complete screening is only an asymptotic law which 
holds in the limit of very large loops with S\p\ — > oo, 
(fT4|l could also be the limiting case of some more general 
diffusion equation which takes into account the contribu- 
tion of other field configurations. In order to derive such 
more general equation, one should make some model- 
dependent assumptions on the form of the cumulants 
of the shifted curvature tensor, so that the coefficients 
j^ai...afc calculated explicitly. In this section 

the diffusion equation ([7|) will be analyzed under the as- 
sumption that the dominating field configurations in the 
vacuum of Yang-Mills theory are thick center vortices, 
which roughly corresponds to the picture of "spaghetti 
vacuum" J^, i2Qi] . Although the assumptions and ap- 
proximations made below by no means capture all prop- 
erties of the theory, they are enough to demonstrate how 
the equation (jl4p emerges in the limit of infinitely thin 
vortices or in the limit of very large loops. Vortices of 
finite thickness can be considered as an effective descrip- 



tion of the interference between thin center vortices which 
were observed in lattice simulations [2l[ and other field 
configurations, as far as such picture of the vacuum of 
Yang-Mills theory correctly reproduces the spectrum of 
string tensions. 

Field strength of a thick center vortex can be obtained 
by smoothing the delta-function in ([TT]) [l8l |: 

i.^) = r {z) e^^^p I dE"/' (y) f{x-y) (21) 
s 

where f {x — y) is some smooth function which decays at 
distances of order of the vortex thickness Ivort- 
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FIG. 2: The probability distribution of L [C; E^] for different 
values of rj (z) S [C] and e [C]: ?7 (z) S [C] = 3, e [C] = 0.05 for 
plot 1, rj (z) S [C] = 1, e [C] = 0.05 for plot 2, r; (z) S [C] = 3, 
e[C] = 0.1 for plot 3 



Proceeding as in the previous section, one can calculate 
the integrals of the cumulants of the shifted curvature 
tensor for the field configurations (|2ip. assuming that 
vortex surfaces E^ are random. However, for thick center 
vortices the linking number L [C; E] can not be rigorously 
defined. Its role is now played by the integral of the 
following form: 

L [C; E] = e^^^p J dS''" (x) J rfE"'' (y) f {x - y) (22) 

S[C] E 

Just as in (|19p . the cumulants of the shifted curvature 
tensor can be expressed in terms of L [C; S]: 



dS'^^-'^ {x,)...dS^^''>'{xk){{F-l,^ (xi). 



pah 



{x,)))= Y.iil'^iC; E,] (z)--- (z) )) (23) 



S[C] S[C] 



On the other hand, the integrals in (|23p can be related to the coefficients r]'^^ - °-k [(7] in ^ using the equation 
These coefficients can be also obtained from P [g; C] using the equation pT|) . Comparing ©, and ITTI) . one can 
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obtain an equation which relates the cumulants of L [C; E^] and the averages of (g): 

S [C] fdgP [g; C] (5) • ■ ■ ^ id) = E (( ^-1 (^) • • ■ T*^' (^) )) (24) 

This relation can be used to find the step distribution P [g; C] if the distribution of L [C; E^] is known. L [C; E^] 
is not in general integer-valued, but it takes integer values when all thick vortices which contribute to the integral 
(|22p are completely inside the loop C. Deviations from integer values occur if some vortices are only partially within 
the loop. If the loop is pierced by right-winding vortices and fc_ left-winding vortices, the average number of 
such vortices can be estimated as Sh [C; E^] ^ {k+ + k-)lyort P [C] /S[C], where P [C] is the perimeter of the loop 
C. For sufficiently large loops and sufficiently small vortex densities one can assume that 6h [C; E^] and k+, fc_ are 
also distributed by Poisson, i.e. that the interaction between vortices can be neglected. As a nearest extension of 
the results obtained in the previous section, one can try to smear the original distribution of linking numbers (|20p 

in such a way that the probability distribution p (h[C; S^]^ is still peaked around integer values, the widths of the 

peaks being proportional to i5L [C; E^]. Without loss of generality the shape of the peaks can be approximated by a 
gaussian distribution. The resulting probability distribution is: 

p (I [C; E,]) = exp {-rj (z) S [C]) ^ ^"^^'^ 1^^^/^^^' ^ - {k+ - k^) , (fc+ + e [C]) (25) 

fe+,fe_ +■ 

where e [C] — Ivort P [C] / S [C\ is assumed to be small and G (x, cr^) = (27rCT^) exp (— x^/ (2fT^)) is the gaussian 
distribution with dispersion a. Such probability distribution for some values of [z) S [C] and e [C] is plotted on Fig. 
[21 It turns out that for the distribution (j25p the generating function for the cumulants can be found explicitly: 



+00 / . -^k 



fc=i 



For comparison, for the distribution (^0)1 W (g, z) ^ rj (z) 5 [C] (cosg — 1). 

For illustration, the step distribution P [g; C] which solves (|24l) will be found for SU (2) gauge group. For SU (2) 
group z = ±1, ^^(1) = and (— 1) = tt^ if the generators Ta are the Pauli matrices. For the cumulants of even 
order the indices ai, . . . , can be again pairwise contracted. As P [g; C] is the function on the group classes, it 
depends only on the absolute value of (g). Taking into account that the Haar measure on SU (2) group classes is 
dg = 2/7rsin^ ^ c?^, one can rewrite the relation (|24p as: 



2/7TS [C] J sin^ e P K; C] iUTrr = (( [C; E_i] )) (27) 


It is also convenient to multiply (P7)) by {—l)'^q'^'' / {2k)l and to sum over all fc, which yields: 

TT 

2/7rS [C] [ d^ sin2 ^ P [C, C] cos {q^/n) = W {q, -1) (28) 



Using the expression (j26p and performing the inverse Fourier transform w.r.t. the variable q — nk, the step distribution 
P [g; C] can be finally found: 

+00 

sin^ ^P[C:C]=-n E cos(fcC)(exp(-7rH2e[C])cos(fc7r)-l) (29) 

k— — oo 



where rj = 77 (—1). over the coordinate ^) looks like a gaussian peak near ^ = 

TT. There is also a (5-function singularity at ^ = 0, which is 
For sufhciently small e (C) the function sin^ ^ P C] irrelevant for diffusion. The width of the peak near ^ = tt 
(which is actually the probability distribution of steps 
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is proportional to \Jlvort P \P\ / S [C]. Thus assuming fi- 
nite thickness of vortices smears the (5-functions in (|14p 
and allows steps by any elements of the gauge group. For 
a fixed loop C the width of the step distribution is pro- 
portional to vortex thickness, but if the vortex thickness 
Ivort is kept fixed, for large loops deviations from equa- 
tion nil) tend to zero as P [C] /S [C] - S [C]"^^^ thus 
asymptotically screening is recovered even for vortices of 
finite size, as it should be [13, HJl • 

The spectrum of string tensions which corresponds to 
the step distribution (j29[) can be found using the equation 
(l9|). As follows from ([9|), string tensions o-r are propor- 
tional to the expectation values of the characters XR id)- 
an = d'^ J dgP [g: C] xr [g)- Direct calculation gives to 
the first order in e [C] : 

(Tk =2r]-2/3rik{k + l) e[C], k = 1/2,3/2,... 

ak = 2/3rik{k + l) e[C], A: = 0,1,2,... (30) 

Remembering that e [C] — IvortP [C] /S [C], one can con- 
clude that the difference between thin and thick vortices 
shows up only in perimeter-dependent effects, which can 
be taken into account if the step distribution P [g; C] is 
allowed to depend on the loop C. Since perimeter depen- 
dence of Wilson loops corresponds to the contribution of 
charge self-energies, it follows from ((50)) that the energy 
required to screen a colour charge in SU (2) representa- 
tion with spin k is proportional to k{k + 1), i.e. to the 
second-order Casimir operator. 

It is interesting to note that Casimir scaling of string 
tensions can also be formally described using the expres- 
sion ([28)) . although the validity of such description should 
be discussed separately. In [24] it was shown that Casimir 
scaling at intermediate distances and asymptotic screen- 
ing can be obtained for Yang-Mills vacuum dominated by 
center vortices if the field distribution (z) / (x) inside 
each center vortex is random at some scale Ic Ivort' 
In this case the distribution of L [C; Sz] should be com- 
pletely different from the distributions (PO)) and ([^5]) if 
the size of the loop is comparable with the vortex thick- 
ness Ivort- By virtue of the Central Limit Theorem it 
should look like a gaussian distribution with the disper- 
sion Sf^l^S [C], where 6f^ is the dispersion of / (x). For 
such distribution W (q) = —Sp l"^ S [C] q^, and the equa- 
tion ([28]) yields the following result: 

sin^ e P [C; C] = -TT^Sf ll fc2 cos (fcC) = 

A;— — oo 

= ^5fll^5{£,) (31) 

The distribution described by (|3ip is singular, but the 
string tensions (7^ can nevertheless be calculated using 




5000 10000 15000 20000 

t 



FIG. 3; Time dependence of 1/2 Tr t; for a particular imple- 
mentation of a random walk on SU (2) group manifold with 
external potential V {g) = —Trg^. 

the equation 

TT 

^k = -^ I sin' C P K; C] Xk iO = 



_ 2d,f ll ^.^^ rf^ sin(2fc + l)$ _ 
^ 2A; -I- 1 gj^ ^ ^ 

^MZ!lfc(fc + i) (32) 

Thus such rather general assumptions concerning the dis- 
tribution of the modified linking number L [C; S^;] lead to 
Casimir scaling at intermediate distances and to screen- 
ing at asymptotically large distances, in full accordance 
with the results obtained in [23 |. 

Finally, it should be noted that as the step distribu- 
tion (|29p is nonzero on the whole group manifold, it still 
describes discontinuous changes of the holonomy g [C] . 
Step distribution which describes continuous changes of 
the holonomy g [C] in nonsingular gauge fields can only 
contain singularities of the form Vai . . . Va^i^ {g, 1) with 
finite k, as in (|3T|) . In order to obtain such step dis- 
tribution for the model of thick center vortices, some 
more advanced assumptions on the distribution of linking 
numbers and on the distribution of fields inside vortices 
should be made. In particular, one could take into ac- 
count that interactions between vortices should lead to 
deviations from Poisson-like distribution of the linking 
number. Although it could be extremely interesting to 
investigate such effects, this is not the aim of this pa- 
per. The derivation presented above is only intended to 
illustrate how the equation (jl4p can arise as the limiting 
case of some more general diffusion equation and how the 
spectrum of string tensions can change in this case. 

VI. CONCLUSIONS 

In this paper a phenomenological analysis of [1, [^, [l3| 

was extended to the case of a general free random walk on 
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the group manifold basing on the exact results obtained 
in [ll| . It was shown that the distribution of steps of a 
random walk of the holonomy g [C] is directly related to 
the spectrum of string tensions of the theory and also to 
certain vacuum expectation values of the curvature ten- 
sor of Yang- Mills fields. In some sense the spectrum of 
string tensions of the theory is a dual image of the spa- 
tial distribution of non-Abelian magnetic flux, just as the 
amplitude of particle scattering in an external potential 
is proportional to the Fourier transform of this potential. 
The description of the dynamics of Wilson loops in terms 
of random walks on the gauge group manifold makes this 
correspondence more apparent. 

The relation between the step distribution of the ran- 
dom walk of g [C] and the spectrum of string tensions 
was used to analyze the random walk of Wilson loops 
in the case of complete screening, when string tension 
an depends only on the A^-ality of the representation R. 
It turned out that when colour charges are completely 
screened, the holonomy g [C] can only change by the el- 
ements of the group center, which implies that the only 
field configurations which contribute to the cumulants of 
the shifted curvature tensor are thin center vortices . 
Complete screening, however, is only an asymptotic prop- 
erty, and the equation (fT4|) can as well be the 5* [C] — > oo 
limit of some more general diffusion equation. This pos- 
sibility was investigated for the case of thick center vor- 
tices in Yang-Mills theory with SU (2) gauge group. It 
was shown that indeed for loops of finite size the equa- 
tion (|14p can only be valid if center vortices are infinitely 
thin, but in the limit S [C] oo it is recovered even 
for vortices of finite thickness. For sufficiently small vor- 
tex sizes or for sufficiently large loops the distribution of 
steps P[g]C] is a gaussian-like peak near g = —I with 
the width e [C] — UortP [C] /S[C]. Of course these con- 
clusions are model-dependent, but at least qualitatively 
the emergence of jumps by the elements of the group 
center can be understood. It could be also interesting to 
study how the equation (jl4p is modified due to contri- 
bution of other field configurations such as instantons or 
monopoles. 

It seems that the most important property of ran- 
dom walks on Lie groups which lead to asymptotical 
screening-like effects is the possibility of fast jump-like 
transitions between points on the group manifold which 
differ by the elements of the group center. For instance, 
if the free diffusion equation ([4]) is modified by introduc- 



ing external Zjv-symmetric potential V (g) , as proposed 
in [1, asymptotically stable probability distribution 
P (g) ~ (ff)) is also ^Ar-symmetric and has N 

equal maxima at some points {gz}, z S Z^. Random 
walker spends most time in the vicinity of such points 
and can only travel from one such point to another in 
a relatively short time. This effect is somewhat simi- 
lar to the restoration of classically broken symmetries 
due to tunneling. In order to illustrate this example, 
random walk on SU (2) group with the external poten- 
tial V (g) = — Tr g^ was simulated. Time dependence of 
1/2Tt g for a particular implementation of such random 
walk is plotted on Fig. [3l It can be seen that indeed the 
random walker mostly wanders near the points g = ±1, 
while the transitions between the regions close to these 
points are very fast and rare. It is reasonable to con- 
jecture that such transitions simulate the discrete jumps 
described by the equation (fT4|) and lead to asymptotic 
screening ^ . 

It could be very interesting to investigate whether such 
jumps could emerge in some explicitly SU (A'^)-symmetric 
dynamical system. It could be also useful to consider 
some generalizations of a free random walk described by 
the equation ([7]). For instance, one could consider ran- 
dom walks with memory, which can be more appropriate 
for the description of domain- like structure of Yang-Mills 
vacuum 0, [11, S HI . 
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